Theory of spin wave excitations of metallic A-type antiferromagnetic manganites 
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The spin dynamics of the metallic A-type antiferromagnetic manganites is studied. An effective 
nearest-neighbour Heisenberg spin wave dispersion is derived from the double exchange model taking 
into account the superexchange interaction between the core spins. The result of inelastic neutron 
scattering experiment on Ndo.45Sro.55Mn03 is qualitatively reproduced. Comparing theory with 
experimental data two main parameters of the model: nearest-neighbour electron transfer amplitude 
and superexchange coupling between the core spins are estimated. 
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Interest in the doped manganese oxides with perovskite 
structure Ri-^B^MnOa (where R is trivalent rare-earth 
and B is divalent alkaline ion, respectively) has revived 
since the discovery of. colossal magnetoresistance (CMR) 
in those compounds. □ 

Extensive experimental study of the doped manganites 
has revealed unusual microscopical structures of spin-, 
charge-, and orbital-ordered states.0 At different dop- 
ing concentration a full range of magnetically ordered 
states such as antiferromagnetic (AFM) insulator, ferro- 
magnetic (FM) metal and charge ordered (CO) insulator 
have been observed. The importance of the charge, spin 
as well as orbitalH and latticetl degrees of freedom and of 
interplay between them has now been well established. 

In the ideal undistorted perovskite structure Mn ions 
are surrounded by six oxygen ions, and form a simple cu- 
bic lattice with the trivalent rare-earth ion sitting in the 
center of the cube. The six oxygen ions, 2 ~, forming 
the octahedra give rise to the crystal field potential that 
removes the degeneracy of d orbital of manganese. As 
a result d level is splitted into t2 g triplet {d xy , d xz , and 
d yz ) and e g doublet (d 3z 2_ r 2, and d x i_ y i). The orbitals 
d 3z 2_ r 2, and d x 2_ y 2 are degenerate and since they point 
toward the O 2- ions they have higher energy then the ti g 
orbitals. The ti g electrons are strongly coupled by intra- 
atomic Hund's coupling in the high spin state S = 3/2, 
and remain localized even in doped compounds. The e g 
electrons are subject to strong Hund's coupling with the 
core spin, they are localized in the undoped compounds 
(one e g electron per manganese ion) due to the strong 
intra-atomic Coulomb repulsion. 

The reference undoped compound LaMn03 is a Mott 
insulator with a spin S = 2, and orbital degrees of free- 
dom. Spins are ordered ferromagnetically in (x, y) plane 
and antiferromagnetically in z direction (so-called A-type 
AFM state). trta The orbitals are also ordered: in-plane di- 
rectional orbitals d 3x 2_ r 2 and d 3l) 2_,.2 alternate in (x,y) 
plane and are stacked in z direction. This leads to the dis- 
tortion of local surrounding of the Jahn-Teller ( JT) active 
Mn ion, rezulting in elongation of MnOg octahedra along 
x (y) direction which has been observed experimentally.Q 

Upon doping of the parent compounds the e g -electrons 



become mobile using oxygen's p-orbitals as a bridge be- 
tween Mn ions and the strong on-site Hund's coupling 
drives core spins to align parallel forming the metallic FM 
ground state. The double exchange (DE) interactions 
that mediates the ferromagnetic coupling between the 
core spins is the key ingredient in a appropriate in- 
terpretation of the phase diagram at the doping range 
0.2 < x < 0.5 where FM metallic behavior is observed. 

Close to the half-doping (x ~ 0.5) the phase diagram 
becomes more subtle. The narrow band manganites (low 
T c compounds) exhibit a charge ordered state close to 
the commensurate filling. This CO state is characterized 
by an alternating Mn 3+ and Mn 4 ions arrangement in 
the (x,y) plane with charge stacking in z-direction. In 
addition, these systems show d 3x 2 _ r 2 / d 3y 2 _ r 2 orbital or- 
dering. In CO state these systems exhibit an insulating 
behavior with a very peculiar form of AFM spin ordering 
at low temperature. The observed magnetic structure is 
a CE-type and consists of quasi one-dimensional ferret 
magnetic zig-zag chains coupled antiferromagnetically.El 
A direct evidence of the CE charge/spin ordered state in 
half-doped manganites has been provided by the electron 
diffraction for Lao.5Cao.5MnO3.EI Similar obs ervat ions 
have also been reported for Pro^Cao^MnOsMIl 1 ! All 
these compounds show paramagnetic insulating states j+t 
high temperature with predominant FM fluctuations.E2l 
Upon lowering the temperature the charge correlations 
together with orbital correlation develops. The FM fluc- 
tuations are first suppressed at charge/orbital ordering 
temperature T Q and completely disappear at the AFM 
ordering temperature Tn .Ej And finally, at T < Tn a CE 
charge/spin/orbital ordered state is established. 

It is clear that all the degrees of freedom involved in the 
formation of ground state are strongly coupled at com- 
mensurate doping. Then the symmetry of the ground 
state is dictated by the competition between the hop- 
ping driven ferromagnetic exchange, antiferromagnetic 
superexchange (SE), electrostatic and lattice elastic en- 
ergies. This competition has been studied in detail to 
interpret the rather complicated CE spin/charge Zrphilal 
ordered state observed in half-doped manganites. Q'ESlJ 

In the moderately narrow band system 
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Ndo.sSro.sMnOalia and Pro. 5 Sro.5Mn0 3 li2 CO state has 
been observed in very narrow band region of Sr concen- 
tration around x ~ 0.5. 

However, in a recent experimental study pf, 
Pro.5Sro. 5 Mn0 3 and Nd .45Sr .55MnO 3 it was foundN 
that instead of the well known CE charge/spin/orbital 
ordered state both systems show A-type AFM order in 
the ground state with uniform d x 2_ y 2 orbital order and 
no clear sign of the charge ordering. The A-type AFM 
order has been also confirmed by the lar.ge anisotropy in 
resistivity and spin wave measurements.EZI 

This finding is consistent with recently obtained 
mean-field phase digram of two-orbital double exchange 
model.Ej It was shown that A-type AFM state and charge 
ordering do not coexist. The A-type spin ordering is un- 
stable against the formation CE-type magnetic structure 
when charge ordering is introduced in the system. In 
the CO state the carrier kinetic energy and hence double 
exchange energy are suppressed, and the CE spin/charge 
ordered state is more favorable due to the additional gain 
in energy by modulating the FM bonds in basal plane and 
generating a "dimerization"-like gap at the Fermi surface 
due to a staggered factor in electron hopping amplitude. 

In the present paper we propose the spin wave theory 
of such an A-type metallic state. We show that the spin 
wave excitation spectrum derived in the leading order of 
1/5 expansion of the double exchange model in A-type 
metallic state reproduces the spin wave dispersion rela- 
tion observed recently in the inelastic neutron scattering 
experiment on Ndo.45Sro.55MnO3.EZl Fitting theory to the 
experimental data gives an estimate of effective transfer 
integral, as well as the AFM superexchange coupling be- 
tween the t2 g core spins. 

In the A-type AFM state electrons are confined in the 
basal plane due to the large Hund's coupling. Since the 
transfer integral between d x i_ y i orbitals is largest in the 
(x, y) plane, the low energy eigen state of two orbital 
tight -binding Hamiltonian has mainly d x 2_ y 2 character. 
Hence there exists an anisotropy in orbital occupation, 
even in the absence of JT coupling. Including JT effect 
will lead to a farther increase of the anisotropy and the 
splitting between d x 2_ y 2 and d^ z 2_ r 2 levels. Therefore, 
in our model we retain only the relevant d x 2_ y 2 orbital, 
and also assume that all the other degrees of freedom are 
integrated out to give effective model parameters. We 
start with the Hamiltonian: 



H = -t [dUjr + Hx 
— Jh ^ Sid; + J SiSj 

m 



fi ^2 Tli 



(1) 



The first term in Eq.(Q) describes an electron hop- 
ping between the the nearest neighbor (NN) Mn-ions in 
(x, y) plane, with t being the transfer amplitude between 
d x 2_ y 2 orbitals. The second term describes the Hund's 
coupling between the spins of localized ti g - electrons Si 



and the itinerant e g electrons with spin oi . The superex- 
change interaction of localized spins between the NN sites 
is given by J, Hi is the particle number operator and fi 
is the chemical potential. 

Next step is to introduce two sublattices A and B with 
spin up and down for the alternating in z-dircction plains 
in A-type AFM state. Then we can expand the core spin 
operators by the Holstein-Primakoff transformation as 
follows: 
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S7 = 



2S_a\ ifieA 
2S0i if i G B 



S — a\ai if i G A 
-S + (3\(3i if i G B 



(2) 



By substituting expressions (g) into the Hamiltonian 
(|l|) and keeping the terms contributing to leading order 
in 1/5 expansion we obtain the following Hamiltonian 
in the momentum space H = H c \ + H sw + H lnt , where 
H e \ and H sw describe the charge and spin degrees of free- 
dom, respectively and interaction among them are given 
by flint: 
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h/3 q /3 q ) + (0 q a q /3- q + H.c. 



where ako- and bk a stands for the electrons on A and B 
sublattice, respectively, £k,t(i) = ^^iWk T JnS/2 — fi 
is the electron dispersion in (x, y) plane including Zee- 
man splitting, 7|| k = l/2(cosfc2; + cosfc y ), t^o, q = Jam + 



7||q), ^q 



J5(l 



-2iq z 



), and m 



2JS - 4J5(1 

l/2(n Qj | — n a .{) = —l/2(nb^ — rib,i) is the magnetization 
of electron subsystem. In the state with fully polarized 
FM planes considered here m = n/2 with n being the 
electron density. 

The spin wave part of the Hamiltonian H svf ([|) is di- 
agonalized by the standard Bogoliubov transformation 



a q = e [u q a q + u q #L q ], /3^ q = 
bringing it into the form 



-Vq&q + Uq/3l q (4) 



a q d q 



(•5) 



The eigen-frequency w q and the coherence factors are 
given by 
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Next we introduce the Fourier transformed retarded 
matrix Green function (GF) for magnons D qw =<C 
A q |A q 3> w where, A q is the two-component operator 

and AJ. = (at, /3- q ). The components of 2? q . w are 
however related among each-other by symmetry relations 
and hereafter we consider the diagonal (normal) (V) and 
non-diagonal (anomalous) (T>) components, expressing 
them in terms of the self-energy operators via the Dyson- 
Beliaev equation 



(7) 



where the following notations has been introduced: 



Dq,u, = \UJ~ A 



•q.ujj 



(8) 
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q,— w 2Sq ?a; 



(9) 



Here E q jW and E q , w are normal and anomalous compo- 
nents of magnons self-energy due to the coupling with 
electrons, -4 q , w is an antisymmetric function in to, while 
and Sq,^ are symmetric functions of u>. 
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FIG. 1. Graphical representation of the normal and anoma- 
lous components of the spin wave self-energy. The black and 
gray straight (wavy) lines stand for the electron (spin wave) 
propagators on A and B sublattices, respectively. 

In the leading order of 1/S expansion, the self-energy 
operators (E and E) are given by the particle-hole bub- 
ble diagrams shown in Fig. 2. The black and gray lines 
stand for the electron propagators on A and B sublat- 
tices, respectively. The upper (lower) graphs corresponds 
to the normal (anomalous) component of the magnons 
self-energy. The anomalous part of the self-energy cou- 
ples degenerate a and (3 magnons denoted by black and 
gray wavy lines, respectively. The dot represents the 
magnon-fermion vertex, analytical form of which is eas- 
ily obtained by reexpressing the interaction term of the 
Hamiltonian (|3j) in terms of transformed magnon oper- 
ators. The above introduced dynamical quantities (|^) 
have the following analytical expression: 
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(10) 



s ± = Jt^K ± V i? Kt - ^k+q|][ek+q| ~ Okf] 



2N 



- [e k+qi - £kt] 2 



where n kcr — [e £kCT + 1] _1 is the Fermi distribution func- 
tion. 

The renormalized spin wave spectrum, arising from 
the second order self-energy corrections, are given by the 
poles of the magnons Green's function (^) , or equivalently 
by the zeros of -D q , w (||). Assuming the large Hund's cou- 
pling Jjj 3> t, and expanding the self-energies up to lead- 
ing order in t/Jjj the spin wave spectrum can be written 
as: 



■ q = 2SJ[J + 2J F (l-7|| q )P- J2 7 2 q , 



(11) 



where 7^ k = cos k z and the effective ferromagnetic ex- 
change energy is given by 

Jf = Jde - J, with J DE = t(a\- l a n )/(2S 2 ) (12) 

being the double-exchange energy given by kinetic energy 
of electrons. 




FIG. 2. The spin wave dispersion of Pro.sSro.sMnOa from 
Ref.17 (the full circles and squares stand for q || (1, 0, 0) and 
q || (0,0, 1), respectively). The solid curves represents fitting 
by Eqjlll including the anisotropy gap. 

The obtained spin wave spectrum ( |Tl| ) corresponds to 
that of nearest-neighbor Heisenberg model with intra- 
plane ferromagnetic exchange integral Jf and antiferro- 
magnetic interplane coupling J. The correspondence be- 
tween spin dynamics in the FM phase of DE model and 
in effective NN Heisenberg model exist at quasi-classical 
level (i.e. leading order in 1/5 expansion) and has been 



studied in Refs. 18 



The long wave-length behavior of the spin-wave dis- 
persion is given by 
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C\\q for q || (1,0,0) 
C ± q for q || (0,0,1) 



(13) 



where Cii = 2SV JJp and C_l = 25 J are intra- and 
inter-plane spin wave stiffness, respectively. Here we 
note, that long wave-length behavior of spin wave spec- 
trum within the FM planes reduces to sound-like disper- 
sion unlike the q 2 behavior inherent for isotropic ferro- 
magnets. 

Fig. 2 shows a fit of inelastic neutron scattering ex- 
periment for Ndo.45Sr .55Mn03 from Ref. [l7| using the 
obtained spin wave spectrum (|ll]) including the phe- 
nomenological anisotropy gap A. The full circles and 
squares stand for the experimental data in (1,0,0) and 
(0,0,1) direction respectively. As has already been 
pointed out in Ref. [IT], the best fit is obtained for a in- 
tralayer FM exchange AJpS ~ 32meV, interlayer AFM 
exchange 2 JS ~ lOmeV, and A ~ 0.26meV. At x = 0.55 
NN correlation function (aj^ajf) = 0.2 and from Eq.[l2] 
one obtains following estimates for d x 2_ y 2 transfer in- 
tegral t ~ 0.195eV and AFM superexchange coupling 
J ~ 3.3meV. 

To summarize, the spin wave dispersion in A-type 
AFM state of the double exchange model is derived tack- 
ing into account the superexchange interaction between 
the core spins. The result of inelastic neutron scattering 
experiment on Ndo.4sSro.55Mn03 is qualitatively repro- 
duced. Comparison of theory with experimental data 
gives the estimates of two main parameters of the model: 
nearest-neighbor electron transfer amplitude t and su- 
perexchange coupling between the core spins J. While 
the present estimate of AFM superexchange strength 
J ~ 3.3meV coincides with that widely used in the 
literaturejia it is factor 4 times larger then SE integral es- 
timated form the Neel temperature-xif the end compound 
CaMnC-3 leading to J ~ 0.8meVE3 In Sr doped com- 
pounds the both e g and t2 g Mn — O transfer amplitudes 
(t e and t g ) are higher than that in Ca doped compounds 
due to the different ionic radii of these element leading to 
different Mn — O distance and Mn — O — Mn bond angle. 
The AFM SE exchange is given by the forth power of the 
transfer amplitude t\ (tg) and hence a small increase of 
the latter might give a larger value of SE strength in Sr- 
doped compounds. This increase of magnetic frustration 
due to the increase of AFM SE strength going from Sr 
to Ca doped compounds has been also observed experi- 
mentallyEil. 
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